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1. Introduction 


The Mittag-LefHer function which is named after the Swedish mathematician 
Mittag-LefHer, who introduced it in 1903 m, has recently been the object of many 
research papers because it plays an important role in Fractional Calculus. It appears 
in solutions of some fractional diHerential and integral equations and as a result it 
has found some physical applications m- 

The Mittag-LefHer function is an entire function of order 1/a, a > 0, defined by 
the power series 


E a (z) 


y _£_ 

h T(i+ka y 


zgC. 


It is easy to see that 


E\{z) = e z , E 2 {z 2 ) = coshz. 

Because the Mittag-LefHer function generalizes the exponential function, one won¬ 
ders if it generalizes or at least shares some properties of the exponential function. 
Several people explored that idea, among them is G. Jumarie, who derived the 
following proposition [8j 
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Proposition 1.1. Suppose that f : M[j~ -A M is continuous and possesses for all 
k £ N Riemann-Liouville fractional derivatives of order ka, with 0 < a < 1. Let 
h > 0. Then 

OO i /ccy 

(1) «*+»)=gfirfh0^ , W’ 

where D a is the modified Riemann-Liouville fractional derivative. 


He used this proposition to show that 

(2) f(x + h) = E a (h a D a ).f(x), 
which he then employed to establish that 

(3) D a E a ( Xx a ) = XE a (Xx a ), x>0. 

This equation is a generalization of the relation 

4-e Xx = Xe Xx , 
ax 

and reduces to it when a = 1. Moreover, he also showed that the semi-group 
property e Xx e Xy = e x ^ x+v ^ may be extended to 

(4) E a (Xx a )E a (Xy a ) = E a (X(x + y) a ), 0 < a < 1, A e C. 

In fact, in jT2] he gave several proofs of this relation. 

Analogous to the definition of cos x and sin x in terms of the exponential function, 
Jumarie defined cos Q (a;“) and sin a (a;“) by 

E a {ix a ) = cos a {x a ) + isin a (x a ), x > 0. 


which leads to 
Therefore, 
is equivalent to 
that is 


cos Q 2 (a; a ) + sin Q 2 (a; a ) = E a (ix a )E a (—ix a ). 
cos a 2 (x a ) + sin Q 2 (x“) = 1, 
E a (ix a )E a (—ix a ) = 1, 


(5) (E^ix*))- 1 =E a (-ix a ). 

But here we note that this relation does not follow from the functional equation 
E a (Xx a )E a (Xy a ) = E a (X(x + y ) a ) 
because if we put y = —x, we have 

E a (Xx a )E a (X(-x) a ) = E a {X{x - x) a ) = E a ( 0) = 1 
which implies that 

(6) (E a (Xx a ))~ 1 = E a (X(-x) a ), 
and when combined with ([5]) leads to 

(7) (^(i® 0 ))" 1 = E a (i(-x) a ) = E a (-ix a ). 
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Under the assumption that there exists a real number M a > 0 such that E a (iM£) 
= 1, the periodicity of the function E a (ix a ), and hence, the periodicity of cos a (a; Q ) 
and sin ct (a"“), will follow since 

E a (ix a ) = E a (ix a )E a (iM%) = E a (i(x + M a )) a ). 


The orthogonality of cos a {mx) a and sin a (m:E)“ now follows from the periodicity 
of these functions together with the fact that 

cos Q 2 (:r a ) + sin Q 2 (:r“) = 1. 

It was also shown that cos a (mx) a and sin Q (ma;)“ form an orthogonal basis for 
L 2 [0,M a ]. In [5], Sec. 6, he obtained an expansion in terms of fractional sine and 
cosine functions as follows. Let / be a periodic function with period M a /u>. Then 
/ can be expanded in a Fourier-type series of the form 

OO 

f(x) = ao/2 + '^2 a n cos a (nujx) a + b n sin Q ,(nwx)“, 

n= 1 


where 

/ \ ol /<jj 

an = 2 \JT ) J f(x)cos a {nujx) a (dx) a , 

and 

/ \ & pM a /cj 

bn = 2 \jr) J f(x)sm a (nu)x) a (dx) a . 

Parseval’s relation takes the form 

f 2 {x){dx) a = i ^{al + bl 

These results paved the way to developing a parallel approach to harmonic anal¬ 
ysis based on these fractional cosine and sine functions. 

Fascinated with these results, the authors of this article tried to derive new results 
using these fractional trigonometric functions and their orthogonality. However, 
they soon realized that there were some errors in Jumarie’s derivations and some 
of the above relations cannot hold. The purpose of this paper is to show exactly 
that. 



2. Preliminary Result 


In this section we show by examples that some of Jumarie’s results are invalid. 
First, we observe from Eq. Q that we should have E a (i(—x) a ) = E a (—ix a ). But 
from the definition 


E a {\x a ) 


\k~otk 
\ -\ A X 

Z^Til + ka ) 1 


we have 


E a (-ix a ) = ^2 

k=0 


(-1 ) k i k x ak 
r(l + ka) 


E a (i(-xD = j2 


(-1 ) ak i k x ak 
r(l + ka) 


k =0 


and 
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and for the last two representations to be equal, we must have (—l) fc = (— l) ak 
which happens only if a = 1. 

We can also show directly that E a (ix a )E a (—ix a ) ^ 1. For, 




£ 

\m —0 


(-ix a ) m 
r(l + am) 


E 


l k,m=0 
oo 


= y iJl ' xan E 


(-1 ) m {ix a ) k+m 
r(l + am)r(l + ak) 

(_l)n-fc 


n—0 
oo 


Ai—0 


r(l + a£;)r(l + a[n — k)) 




n —0 


where 


A n(a) = Y 


(- 1 ) 


i—k 


k—0 


r(l + afc)r(l + a(n — k))' 


It is easy to see that A 0 = 1, A 2k+ 1 = 0, Vfc e N 0 . However, 

1 1 1 

2 ~~ r(l + 2 a) r 2 (l + a) + T(l + 2 a) 
2 1 

r(l + 2a) ~ r 2 (l + a) 

2r 2 (l + a) — r(l + 2a) 
r(l + 2 a)r 2 (l + a) ^ 


In fact, all the even A 2k 7 ^ 0, but all of them are equal to zero only if a = 1. 

In Figure[L] the product p(x, a) := E a {ix a )E a (—ix a ) is plotted for several values 
of 0 < a < 1. 



Figure 1. The product E a (ix a )E a (—ix a ) for several values of 
0 < a < 1. 
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Moreover, if there existed a real number M > 0 such that E a (iM a ) = 1, we 
must have cos„(M“) = 1, sin Q (M“) = 0. But for example, for a = we have 

00 T k+ 1/2 

cos 1/2 (a; 1/2 ) = e~ x and sin 1/2 (a; 1/2 ) = ^(~ 1 ) fc , v 

Using the Pochammer notation, we have T(k + 3/2) = T(3/2) (3/2)& and hence 
sini /2 (a; 1/2 ) = 


where $( 0 , 6 ; x) is the confluent hypergeometric function. Using the relation 

<l>(a, b; x) = e x $(6 — a, b; — x) 

we obtain 

$(1,3/2; -x) = e _x $(l/2,3/2; a;), 
or 

r l/2 

sini / 2 (a: 1/2 ) = ^-^ye“ a: $(l/ 2 ,3/2; x). 

But in view of the relationship between the confluent hypergeometric and the error 
function 

c z 

Erf(z) = 0 $ (1/2,3/2; -z 2 ) = / e^dt, 

Jo 

it is clear that the only solution for cosi/ 2 (Af 1 / 2 ) = 1 and sin 1 / 2 (M 1 / 2 ) = 0 is 
M = 0, i.e., there is no period for these fractional trigonometric functions for 
a 2 - 

In the next section we show that our results hold not only for a = 1/2 but for 
all 0 < a < 1 . 


r 1/2 


r (3/2) 


HUD 1 


k—0 


(3/2)* 


r l/2 


E 


r (3/2) ^ (3/2)*;(l)fc 

x 1/2 y. (l) k (-x) k _ 
B(3/2)^ (3/2 ) k k\ ~ T(3/2) 


r l/2 


$(1,3/2;-*), 


3. The Main Result 

In this section we extend the results of the previous section from a = 1/2 to all 
0 < a < 1. To this end, we need the following definition. 

Definition 3.1. A function f : R + — > R is called completely monotonic (CM) if 
f eC°° and Vn € N 0 Vx € R+ : (-1 )"/W(*) > 0. 

Theorem 3.1. Let 0 < a < 1. Then there does not exist an M a > 0 so that 
E a (iM(f) = 1 or, equivalently, cos a (M") = 1 and sin a (M“) = 0. 
Moreover, the relation 

(8) E a (\x a )E a {\y a ) = E a (X(x + y) a ), 0 < a < 1, A £ C, 

can hold only if a = 1. 
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Proof. We first note that cos a (x a ) = E 2 a {—x 2a ). This is verified by direct compu¬ 
tation or by referring to the duplication formula m Eqn. 2.14]. It suffices to show 
that cos a (x“) < 1 for all x > 0. 

To this end, we use a result from [IS] which states that the Mittag-LefHer function 
E a {—x) is completely monotonic for all 0 < a < 1 and x > 0. As x H► x a is a 
Bernstein function for 0 < a < 1 and the composition E a {— ■) o (•)“ of the CM 
function E a (— •) with the Bernstein function x H > x a is again a CM function m 
Theorem 3.7. (ii)], we see that E 2 a (—x 2a ) is a CM function for 0 < a < 

As E 2a ( 0) = 1 and E 2a {—x 2a ) —> 0+ (see, for instance, [TSJ Eq. (3.3)]), the 
complete monotonicity for 0 < a < 1 and x > 0 now implies that there cannot 
exist an M a > 0 such that cos Q (M“) = 1. 

Now suppose | < a < 1. It is shown in m Section 5.2] that in this case the 
Mittag-LefHer function E 201 can be written as a sum of two functions, / 2 a and g 2 a , 
where the former is CM and the latter oscillatory. More precisely, 

E 2a {~x 2a ) = f 2 a ( x 2a ) + g 2 a {-x 2a ) 1 i < a < 1, x > 0, 

where 

2q . _ 1 [°° e~ sx2a s 2a sin(2a7r) 

2a x 7i Jo s 4a + 2s 2 “ cos(2a7r) + 1 

and 

32 «(-z 2 “) := | e~ x2a cos (-/2“) cos (x 2 “ sin . 

The function f 2a is CM, satisfies f 2a (0) = 1 — and increases towards zero from 
above, whereas the function g 2a is oscillatory with exponentially decaying amplitude 
and <? 2 a(b) = Hence, E 2a {—x 2a ) < 1, for all | < a < 1 and x > 0. Thus, there 
cannot exist an M a in this case either. 

The case a = \ was considered above. 

Next we show that the purported functional equation (jSJ can hold only if a = 1. 
First note that E a (Xx a ) = E a o h(x), where h(x) := Xx a . Writing F a for the 
composition E a o h, Eq. © reads 

(9) F a (x)F a (y) = F a (x + y). 

However, the only nonzero continuous solutions of Eq. (|UJ| are exponential functions 
of the form e cx , where c £ C; see for instance [TJ Chapter 2], Therefore, F a (x) = 
E a (Xx a ) = e cx , for all x. Successively differentiating the power series for E a (Xx a ) 
and e cx and letting x —> 0, shows that A = c and a = 1. □ 


Figure [2] below depicts some graphs of cos Q (x“) = E 2a {—x 2a ), for 0 < x < 1. 
The case a := 1 produces the cosine function. 

In [T2] . the purported functional equation for E a (Xx a ) is derived using several 
approaches. One of them is the product rule 

(10) D a (fg)=gD a f + fD a g. 


However, this rule is not correct. 

To see this, let g := / := x 1 / 2 . It is straight-forward to establish that for x > 0, 
p > — 1, and 0 < a < 1, the fractional derivative of x p is given by 


D a x p = 


r(i + P ) 
r(l +p — a) 
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Hence, the left-hand side of m computes to 


D a {x) 


x 1 “ 

r(2 -ay 


whereas the right-hand side equals 

Both sides are identical only if a = 1. 

Similarly, the two purported chain rules, which are also employed in the deriva¬ 
tion of ©, namely 


D a (/ o u) ( x ) = D a u(x ) 

du 


and D a (fou)(x) = (DZf)-(^ 


are not correct. To validate this, take for the former, f{x) := x 2 and u{x) := x 1 / 2 , 
x > 0. Then, for any 0 < a < 1, one has 

(11) D a (f ou){x) = D«x= f " 

T(2 — a) 

and 


( 12 ) 


df 


du 


U=X 1 / 2 


D a u(x) 


2x 1/2 D a x 1/2 


y/nx 1 “ 

r (!-a)’ 


and this two expressions are identical only if a = 1. Now for the latter, take 
f(u) := u x ! 2 and u(x) := x 2 , x > 0. Then, for any 0 < a < 1, we again have 


D a {fo U )(x) = D a x = 


whereas now 


T(2 — a) ’ 

(is) ( du ) u, ■ (|)“ = r>;<» V2 )L., a (2x)“ = 


2°‘- 1 ^x 1 - a 

r(|-a) 


Eqns. ©, and m are identical only if a = 1. 

For illustrative purposes, we plotted the difference E a {\(x+y) a )—E a {\x a )E a {\y a ) 
in Figure [5] for (x,y) € [0,2] 2 , A = 1, a = j (left), and a = | (right). 
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Figure 3. The difference E a {X{x + y) a ) — E a (Xx a )E a (Xy a ). 
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